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Figure 1 : Lowest order contribution to Bs — Bg mixing in the Standard Model. 



1 Introduction 

Currently the prime focus of experimental elementary particle physics is the investi- 
gation of the flavor sector of the Standard Model. Transitions between different fermion 
generations originate from the Higgs- Yukawa sector, which is poorly tested so far. The ex- 
perimental effort is not only devoted to a precise determination of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix which parameterizes the flavor-changing couplings. Flavor- 
changing neutral currents (FCNC) also provide an ideal testing ground to search for new 
physics, because they are highly suppressed in the Standard Model: FCNC's are loop- 
induced, involve the weak coupling constant and the heavy W boson, are suppressed by 
small CKM elements or the GIM mechanism and further often suffer from a helicity- 
suppression, because flavor-changing couplings only involve left-handed fields. Therefore 
experiments in flavor physics are much more sensitive to new physics than the precision 
tests of the gauge sector performed in the LEP/SLD/Fermilab-Run-I era. Decays of B 
mesons are especially interesting: they allow us to determine three of the four CKM pa- 
rameters, their rich decay spectrum helps to overconstrain the CKM matrix, they have the- 
oretically clean CP asymmetries (as opposed to K irn decays), information from Bd, Bs 
and 5+ decays can be combined using SU (3)i? symmetry, the large b quark mass permits 
the use of heavy quark symmetries and the heavy quark expansion, and in many extensions 
of the Standard Model third generation fermions are most sensitive to new physics. 

While Bs mesons cannot be studied at the B factories running on the T(45') resonance 
[^, they are copiously produced at hadron colliders Bg mesons mix with their antipar- 
ticles. Therefore the two mass eigenstates Bh and Bl (for "heavy" and "light"), which are 
linear combinations of Bs and Bs, differ in their mass and width. In the Standard Model 
Bs — Bg mixing is described in the lowest order by the box diagrams depicted in Fig. |T]. 
The dispersive part of the Bs — Bs mixing amplitude is called M12. In the Standard Model 
it is dominated by box diagrams with internal top quarks. The absorptive part is denoted 
by and mainly stems from box diagrams with light charm quarks. is generated 
by decays into final states which are common to Bg and Bg. While M12 can receive siz- 
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able contributions from new physics, is induced by the CKM-favored tree-level decay 
b CCS and is insensitive to new physics. Experimentally Bg — Bg mixing manifests it- 
self in damped oscillations between the Bg and Bg states. We denote the mass and width 
differences between Bh and Bihy 



Am = Mh- Ml 



AT 



By solving the eigenvalue problem of M12 — 1^12/ 2 one can relate Am and AF to M12 and 
ri2: 



Am = 2 |M- 



12 



Ar = 2 |ri2| cos< 



where (p is defined as 
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(2) 



Am equals the Bg — Bg oscillation frequency and has not been measured yet. In deriving 
(|l[) terms of order |ri2/Afi2p have been neglected. in (Q) is a CP-violating phase, which 
is tiny in the Standard Model, so that AFsm = 2|ri2|. Unlike in the case of B^ mesons, the 
Standard Model predicts a sizable width difference AF in the Bg system, roughly between 
5 and 30% of the average total width T = (F^ + Th)/2. The decay of an untagged Bg 
meson into the final state / is in general governed by two exponentials: 



F[/,t]oce-r-*|(/|5L 



f + e-^"'\{f\BH[ 



(3) 



If / is a flavor-specific final state like -D^tt^ or Xi^u, the coefficients of the two exponen- 
tials in (|[) are equal. A fit of the corresponding decay distribution to a single exponential 
then determines the average width F up to corrections of order (AF)^/F. In the Standard 
Model CP violation in Bg—Bg mixing is negligible, so that we can simultaneously choose 
Bi and B^ to be CP eigenstates and the b ccs decay to conserve CP. Then Bh is CP-odd 
and cannot decay into a CP-even double-charm final state fcp+ like {J/i'(p)L=o,2, where L 
denotes the quantum number of the orbital angular momentum. Thus a measurement of the 
Bg width in Bg fcp+ determines Tl- By comparing the two measurements one finds 
AF/2. CDF will perform this measurement with Bg Djn^ and Bg J/i'(p in Run-II 
of the Tevatron . 



2 QCD corrections 

Weak decays of B mesons involve a large range of different mass scales: first there is 
the W boson mass Mw, which appears in the weak b ccs decay amplitude. The second 
scale in the problem is the mass m^, of the decaying b quark. Finally there is the QCD scale 
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parameter Aqcd, which sets the scale for the strong binding forces in the Bg meson. QCD 
corrections associated with these scales must be treated in different ways. To this end one 
employs a series of operator product expansions, which factorize the studied amplitude into 
short-distance Wilson coefficient and matrix elements of local operators, which comprise 
the long-distance physics. Here in the first step the W^-mediated b ccs decay amplitude 
is matched to matrix elements of local four-quark operators. We need the two | Ai?| = 1 
current-current operators 

Qi = Q7m(1 - l5)bjSjYi'^ - lb)ci Q2 = Ci7^(l - 75)6jSj7^(l - 75)0^-, (4) 

where i,j are color indices. Q2 is pictorially obtained by contracting the W line in the 
b CCS amplitude to a point. Qi emerges, once gluon exchange between the two quark 
lines is included. In the effective hamiltonian 



Gf ^ 



(5) 



the Wilson coefficients Cr are determined in such a way that the Standard Model amplitude 
is reproduced by ( ccs iHefj] b) up to terms of order m^/M^. The Fermi constant Gp 
and the CKM elements have been factored out in @. The C^'s contain the short-distance 
physics associated with the scale Mw. QCD corrections to the Wilson coefficients can be 
computed in perturbation theory. The renormalization group evolution of the C^'s down to 
the scale fii = 0{mh) sums the large logarithms \n{iii/M]v) to all orders in perturbation 
theory. The minimal way to do this is the leading log approximation which reproduces all 
term of order a" ln"(yUi/Mvi/), n = 0, 1, . . ., of the full Standard Model transition ampli- 
tude. The next-to-leading order (NLO) corrections to the coefficients comprise the terms 
of order a""*"^ ln"(/ii/Miy) and have been calculated in We remark that there are also 
penguin operators in the effective hamiltonian Tig//- We have omitted them in (Q), because 
their coefficients are very small. Their impact is discussed in 
ArsM = 2|ri2| is related to Heff by the optical theorem: 

1 
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(6) 



Here 'Abs' denotes the absorptive part of the amplitude, which is obtained by retaining 
only the imaginary part of the loop integration. The corresponding leading-order diagrams 
are shown in Fig. ^ In the next step of our calculation we perform an operator product 
expansion of the RHS in @ in order to describe in terms of matrix elements of local 
\AB\ = 2 operators: 



|Abs {Bs\i J d^xTHeff{x)rLeff{Q)\Bs 



Gp^b |T/*T/ |2 

F ( 41 {Bs\Q\Bs) + Fs('^] WQslBs) 



mt 
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QCD 

rrib 



(7) 
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Figure 2: Leading-order diagrams for F 
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The two dimension-6 operators appearing in are 

Q = Si7/^(1 - l5)biSjY{l - j5)bj, Qs = Si{l + 75)6iSj(l + 75)6^- 



(8) 



In the leading order of QCD the RHS of (0) is pictorially obtained by simply shrinking the 
(c, c) loop in Fig. ^ to a point. Our second operator product expansion is also called heavy 
quark expansion (HQE), which has been developed long ago by Shifman and Voloshin [^. 
The new Wilson coefficients F and Fs also depend on the charm quark mass rric, which is 
formally treated as a hard scale of order irib, since rric ^ Aqcd- Strictly speaking, the HQE 
in (0) is an expansion in Kqcd/ \Jml — 4m^. For the calculation of F and Fs it is crucial 
that these coefficients do not depend on the infrared structure of the process. In particular 
they are independent of the QCD binding forces in the external Bg and Bg states in so 
that they can be calculated in perturbation theory at the parton level. The non-perturbative 
long-distance QCD effects completely reside in the hadronic matrix elements of Q and Qs- 
It is customary to parametrize these matrix as 



{Bg\Q{^^2)\Bs) 
{Bg\Qs{ii2)\Bs) 



-JlMlB{iX2) 



Ml 



- g fl Ml - — - — — — Bs{ii2) 
3 [mb{^2) + ms{^i2)Y 



(9) 



Here Mb^ and fs^ are mass and decay constant of the Bg meson. The quark masses m?, and 
nig in (|^) are defined in the MS scheme. In the so called vacuum insertion approximation 
B(fi2) and Bs{n2) are equal to 1. /i2 = 0{mh) is the scale at which the \/S.B\ = 2 
operators are renormalized. It can be chosen different from fii. The dependence of AF 
on the unphysical scales fii and /i2 diminishes order-by-order in perturbation theory. The 
residual dependence is usually used as an estimate of the theoretical uncertainty. The /ii- 
dependence cancels between the | Ai?| = 1 Wilson coefficients Ci 2 in @ and the radiative 
corrections to F and F5 in (^. The terms in F and Fs which depend on fi2 cancel with 
corresponding terms in -B(/i2) and -Bs(/^2)- The scale /i2 enters a lattice calculation of these 
non-perturbative parameters when the lattice quantities are matched to the continuum. 

The leading-order calculation of AF requires the calculation of the diagrams in Fig. || 
and has been performed long ago [|T0|]. Subsequently corrections of order Agco/'mb to 
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Figure 3: Next-to-leading-order diagrams for ri2. Qs is the chromomagnetic penguin operator. 

have been computed in The next-to-leading order calculation requires the calculation 
of the diagrams depicted in Fig. ^ The motivations for this cumbersome calculation are 

1) to verify the infrared safety of F and Fs, 

2) to allow for an experimental test of the HQE, 

3) a meaningful use of lattice results for hadronic matrix elements, 

4) a consistent use of A^, 

5) to reduce the sizable -dependence of the LO, 

6) the large size of QCD corrections, typically of order 30%. 

The disappearence of infrared effects from the Wilson coefficients F and Fs mentioned 
in point 1) is necessary for any meaningful operator product expansion. Yet early critics 
of the HQE had found power-like infrared divergences in individual cuts of diagrams of 
Fig. |[ In response the cancellation of these divergences has been shown [[11]], long ago 
before we have performed the full NLO calculation. However, there are also logarithmic 
infrared divergences. We found IR-singularities to cancel via two mechanisms: 

• Bloch-Nordsiek cancellations among different cuts of the same diagram. 
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• factorization of IR-singularities, which end up in {Bs\Q\Bs), {Bs\Qs\Bs) ■ 

Point 2) above addresses the conceptual basis of the HQE, which is sometimes termed 
quark-hadron duality. It is not clear, whether the HQE reproduces all QCD effects com- 
pletely. Exponential terms like exp(— Kmfe/AQCD)^ for example, cannot be reproduced by 



a power series [|T2p. The relevance of such terms can at present only be addresses exper- 
imentally, by confronting HQE-based predictions with data. The only QCD information 
contained in the LO prediction for AF is the coefficients of a" In" Mw, associated with 
hard gluon exchange along the W^-mediated h ccs amplitude. The question of quark- 
hadron duality, however, addresses the non-logarithmic QCD corrections, which belong to 
the NLO. While it is certainly very interesting to find violations of quark-hadron duality in 
B physics, it will be hard to detect them in AF: in the LO diagrams in Fig. ^ the heavy c,c 
quarks recoil back-to-back against each other and are fast in the b rest frame. The inclusive 
b CCS decay is more sensitive to uncontrolled long-distance effects, because in some 
parts of the phase space the c,c quarks move slowly in the b rest frame or with respect to 
each other. Still the HQE prediction for b ccs [[T3p agrees with experiment [[T^]. If 
one further takes into account that AF has an overall hadronic uncertainty associated with 
fs^B and fs^Bs, it appears very unlikely that violations of quark-hadron duality can be 
detected in AF. Points 3) and 4) are related to the fact that leading-order predictions are not 
sensitive to the renormalization scheme, which impedes the lattice-continuum matching of 
the non-perturbative parameters. Likewise the /X2 dependence of this matching procedure 
cannot be addressed in the leading-order. The yUi-dependence of AF is huge in the leading 
order. It is reduced in the NLO, but still remains sizable. The results for F and Fs can be 
found in Tab. ^. The reduction of the fii dependence can be verified from the table. The 
numerical values of the NLO coefficients depend on the renormalization scheme. The pre- 
cise definition of this scheme involves the subtraction prescription for the ultraviolet poles 



(dimensional regularization with MS subtraction [[T5|]), the treatment of 75 (for which we 
have used the NDR scheme) and the chosen definitions of the evanescent operators [[T^, 
which can be found in [|]]. The lattice-continuum matching must be done in the same 
renormalization scheme, so that all scheme dependences cancel in the prediction for AF. 
Including the corrections of order AQco/mb [0] our NLO prediction reads 



AFsM ^ / fB. y 
F I 245 MeV / 



(0.234 ± 0.035) Bsinib) - 0.080 ± 0.020 ] . (10) 



Here mb{mb) + ms(mb) = 4.3 GeV (in the MS scheme) and ml/nil = 0.085 has been 
used. Since F is small, the uncertainty in B is irrelevant, and the term involving FB has 
been absorbed into the constant -0.080 ± 0.020 in Recently the KEK-Hiroshima 
group succeeded in calculating /b^ in an unquenched lattice QCD calculation with two 
dynamical fermions [[T7|]. The result is fs^ = (245 ± 30) MeV. A recent quenched lattice 
calculation has found Bs{mb) = 0.87 ± 0.09 [ P^ ] for the MS scheme. A similar result has 
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Table 1: 





mb/2 


rrib 


2mfe 


-Fs 


0.867 


1.045 


1.111 




1.729 


1.513 


1.341 


F 


0.042 


0.045 


0.049 


i?(o) 


0.030 


0.057 


0.103 



Table 2: Numerical values of the Wilson coefficients F and Fs for ml /ml = 0.085. Leading- 
order results are indicated with the superscript (0). The precise definition of our renormalization 
scheme can be found in llsll. 



been obtained in [[T9|]. With these numbers one finds from ([T0|): 

= 0.12 ±0.06. (11) 
Here we have conservatively added the errors from the two lattice quantities linearly. 



3 New physics 

In the presence of new physics arg M12 and thereby (p in can assume any value. 
Non-standard contributions to (p can be measured from CP-asymmetries, which requires 
the resolution of the rapid Bg—Bs oscillations and tagging, i.e. the discrimination between 
Bs and Bg mesons at the time t = of their production. From (|l[) one verifies that a 
non-vanishing also affects AF, which can be measured from untagged data samples and 
therefore involves better efficiencies than tagged studies. Of course in the search for new 
physics AF is only competitive with CP asymmetries, which determine sin 0, if is not 
too close to or ±7r. Nevertheless the information on from both tagged and untagged 
data should be combined. 

As discussed at the end of Sect. [I], AF is most easily found from the lifetimes measured 
in the decays of an untagged Bg sample into a flavor-specific final state and into a CP- 
specific final state fcp, respectively. In the presence of a non-zero CP-violating phase 
the mass eigenstates B^ and Bh are no more CP eigenstates, so that now both exponentials 
in (Q) contribute to the decay Bg fcp. Then this method determines 

AF COS0 = AFsM cos^ 0. (12) 

As first pointed out in [pO|], one can determine | cos0| without using the theoretical input 
in ([T0|): if one is able to resolve both exponentials of (|^) in the time evolution of a Bg 
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decay into a flavor- specific final state, one will measure the true |Ar|. By comparing 
with ( [T^ ) one can then solve for | cos0|. This method, however, requires to distinguish 
cosh((Ar)t/2) from 1 and is very difficult to carry out. In ^A\\ a different method has 
been proposed, which only requires to measure lifetimes and branching ratios: first define 
CP eigenstates Bf^ and 5f such that Bf'^ ^ DtDj. Then define 



AFcp is related to as 



ATcp = r(5r")-r(fin. (13) 



AFcP = 2|ri2| 



Hence AFcp equals AFsm, but is not affected by the new physics phase (j) ^U! By 
measuring both AFcp and AT cos one can infer | cos0| from (|T2|). Loosely speaking, 
AFcp is measured by counting the CP-even and CP-odd double-charm final states in Bg 
decays: 



, (-) 



AFcp = 2F ^ Br { B f) (l - 2xf) 



(14) 



Here Br ( / ) is the branching ratio of an untagged Bg meson into the final state /, 

F is the average Bg width, the sum runs over all double-charm final states and xj is the 
CP-odd component of the final state /, e.g. xj is for a CP-even state and equals 1 for a 



CP-odd state. In the Shifman-Voloshin limit one can show that AFcp is exhausted by 



the D^*^'^D^*^ final states [^. Moreover these four final states are purely CP-even in this 
limit. ALEPH has measured the sum of these branching ratios [ |24l ] and found, relying on 
the SV limit, 

AFcp ~ 2 5r ( dW+^W" ) = 0.2Qtlf,. (15) 

In the future one can extend this method by including all detected double-charm final states 
into the sum in ( [14] ) and determine the CP-odd fraction Xf of each final state by measuring 
the Bg lifetime in the studied mode [p]]. 
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